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ABSTRACT 

A computer program has been developed for determining the 
subsonic pressure, force and moment coefficients for a fuselage- 
type body using slender body theory. The program is suitable 
for determining the angle of attack and sideslipping character- 
istics of such bodies in the linear range where viscous effects 
are not predominant. Procedures have been developed which are 
capable of treating cross sections with corners or regions of 
large curvature. 
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SUMMARY 


A computer program has been developed to determine the subsonic 
pressure^ force and moment coefficients based on slender body theory for 
bodies of arbitrary cross section^ The program is based on the integral 
repiesentation of the potential in which the flow in the crossflow plane io* 
considered to be induced sources distributed about the cross sectional 
boundary. Analytical expressions are derived for cd and its derivatives 
and the integrals appearing in these are evaluated by dividing each cross 
sectional boundary into straight line segments approximating the integro.nds 
over these segments. Results for pressure force and moment coefficients 
have been obtained for circular cone and ogive bodies and compared with 
analytical determinations from slender body theory. Results arc also 
obtained for a typical slab-sided** fuselage. 

In Part III modifications have been developed which extend tho 
applicability of the program in Pa t II to crossections with corners or local 
regions of high curvature. 
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INTRODUCTION 


Computerization of aerodynamic theory has progressed to a point 
where the flow field analysis of complete aircraft configurations by a 
single program is now an attainable goal* Programs designed for*£his 
purpose do in fact exist, but predictably they are extremely large and 
abound with subtleties often not evident to the user. Generally, each new 
application undergoes a '^debugging*' stage which may in iu elf constitute 
a major effort. Much of the complexity of these programs is attributable 
to the level of precision of the underlying theory. Although often impres- 
sive, this precision is not always required. In some stages of preliminary 
design, for instance, it would be more desirable to sacrifice precision for 
simplicity. One approach in this spirit is to replace the commonly 
employed exact superposition method which panels the entire aircraft sur- 
face, placing appropriate singularities at each panel, with linearized 
theories involving only solutions of a local two-dimensional potential equa- 
tion. In the exact theories a determination of the singularity strengths 
required to satisfy boundary conditions leads to the necessity of inverting 
very large matrices. The quasi-two-dimensional nature of linearized 
theories on the other hand considerably reduces the size of the matrices 
encountered and consequently places far less demand upon computer capa- 
bilities. 

It is the purpose here to develop programs based ot slender body 
theory, utilizing two-dimensional singularities distributed along a cross 
sectional contour to solve for the required potential function in the cross 
flow plane. Such an approach is felt to be particularly adaptable to the 
formulation of the interaction problems encountered in the analysis of 
complete aircraft configurations. 
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SYMBOLS 



g(x) 


h 

i 

i, i+4l 


iL 




M 



n(i, n) 


N 


Coefficient of doublet term in expansion of complex Potential W. 
Cross sectional boundary at station x . 

Lift and side force coefficients. 

Pitch and Yaw moment ■coefficients aboTzt nose -of body. 

Cross sectional boundary at x = . 

Pressure coefficient (p - p^)/(pU®/2) 

Horizontal and Vertical Force components for body of unit length. 
Function of x derived from outer solution to potential equation. 
Radius of curvature of cross sectional boundary. 

Index of points along cross sectional boundary C. 

Segment of C from i to i+1. 

Total number of segments into which C is divided. 

Length of segment 4, i+1 on C^. 

Mach number 

Components of moments about nose for body of unit length. 

Inner unit normal to segment i,i+l.- 
Total number of stations 
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p 

Po 

q* 

R(i, j, n) 
lR(i. j. n) 
s 
S 

u(i, n) 

U 


W 

j£ 

y»2 

z 

OL 


Pressure 

Stream Pressure 

V® + w® 

Displacement from pt P. ^ to pt P/ ^ on C^. 

Vector displacement from P. to P/ 

^ j,n 

Distance along C^. 

Cross sectional area. 

Unit tangent to segment i,ii-l. 

Fiee Stream Velocity. 

Normalized Radial Polar Coordinate. 

Normalized y component of velocity in wind azes. 
Normalized z component of velocity in wind axes. 
Normalized Complex Potential Function. 
Normalized Longitudinal Function. 

Wind axe if coordinates in transverse plane, 
y + iz 

Complex location of cross sectional centroid. 
Angle of attack . 
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0 

6 ( ) 

e 

a 

K 

6 

T1 

cp 

$ 

6v /6x 
o 

6v/fix 

c 

¥ 


Vl - 

Differential corresponding to displacement normal to C^. 

Angle subtended by a, i+L at pt. j at station . (see Fig. 3) 

Angular Polar coordinate. 

2 Dimensional source density. 

Value of Q at point P. 

^ 1 , n 

Angle between tangent to C and y axis 

Normal displacement from mid point of i, i+1 on to i, i+1 

Pex'turbation potential, 
cp i g(x) 

Ucp^ 

Body slope in body axis frame of reference. 

Body slope in wind aucis frame of reference. 

Complex position on C in wind axes. 

Complex position on C in body axes. 

Yaw angle. 
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PART I 


THEORY AND DEVELOPMENT OF NUMERICAL PROCEDURES 

A. Syi.opsis of Subsonic Slender Body Theory 

According to slender bcdy theory (ref. 1), the flow dtsttl^banc6L^ • - 
near a stifficiently regular 3-D body may be represented by a potential 
in the form: 

f(xyz) = Ucp^ = U[cp(xyz) + g(x)] (1) 

cp(xyz) is a solution of the 2-D Laplace equation in the y, z cross flow plane 
satisfying the following boundarv conditions appropriate to wind axes* 

^ = 0 at ® (2a) 

C(x), n, and being defined in Fig* (1)* A general solution for Cp may 
be written as the real part of a complex potential function W(Z) with 
Z = y + iz. 

CD = ReW = RcFa (x) InZ + Y A (x)/z"l (3) 

* e o n -J 

n 

A useful alternative representation of co and W is obtainable with the aid 
of Gri|)n's theorem, (ref. Z) 

cp=ReW=-2Re (f 0(0 In (Z-C)ds (4) 

c(x) 

where 0(C) is a "source' density for values of C • y + iz , (y , ■ ) 

c c c c 

being coordinates of a point on the contour c(x). 


Although windJ^es have been adopted as a reference, the computations 
have been forintxlated in terms of input data obtained from a body axes 
frame of reference. This avoids the necessity of generating new input 
data for each change in bouy c i'ltude. 
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The function g(x) is obtained by matching $ of Eq, (1) which is 
%alid in the neighborhood of the body with an appropriate "outer' solution, 
gfx) is then found to depend explicitly on the longitudinal variation of 
cross sectional areas S(x) , i.e.: 


g(x) 


f X 1 

= ^ S'(x)ln(P/2) - ^ J S'(t)ln<x-t)dt+l Ji 

L c X 


S"(t) ln(t-x)dt 


S'(0) 


Inx 


1 


S(l) 

2 


In(l-x) 


^ = yi.M® 


(5) 


The pressure coefficient, to an approximation consistent v ith 
slender body theory is given by the expression: 



The force arsd moment about the origin on the portion of the body between 
the nose and station x are represented by the coefficients: 


5= 2iJA.(x) +;^(s(x) Z (x)j 

PU* 1 dx \ g / 




M ♦iM 

y 2 

pu^ 


X 


. F +ir 

i^.JL 5 _ 2n I A.(t)dt - S(x)Z (x>|- 

f pTj2 •' ! g ^ 


<3) 


where Z (x) = T + iz reoresents the complex location of the cross 

g » g * 

sectional centroid at station x, and A^(x) is the coefficient of th.' 1/Z 


term of (S), In terms of these force and moment expressions the 


Of i'uui. iH'AU'nr 
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HM>t« commonly used aerodynamic coefficients are written; 



.Ms ,3 

c .. = K-^)t 


ref ®fef 


,M s 


- ^ s . 

ref 


where I =. body length and ^ref convenient reference length and 

area respectively, usually, determined by the overall configuration tc be 

analyzed. For this report has been chosen to be equal to L and 

S f=jJ. 

ref 


The reduction of cofr.putnfions of these expressions to a numerical 
proceduiC shall be based on the integral representation of given in Eq. (4). 
The point of departure sha’l be the discretization of le cross sectional 
bmsndary into a large number of short linear segments over each of which 
the source density 0 shall be assun^^d conjtant at a value to be determined 
by houtidary cofiditions. 



fiw Sommary of Equations. Compatational Procedures and Sample 

^ 

Derivations of the equations presented in this section are given in 
Appendix A. 

Since analytical results for bodies of revolution are readily available 
computations have been carried out for the purpose of comparison in the 
cases of a circular cone; 

r(x) 3 X tan 10® . 0 < x < 1 

and an "ogive" of circular cross section; 

r(x) = x(l - |) tan 10® 0 <x < 1 


n 


both at angle of attack a = . 1 and at zero Mach no. 

1. Processing of Surface Data 

The original data consists of the cross s'ectional boundaries C_ 

at each presented in body axes coordinates as shown in Fig. 2. 

Starting at a convenient station x^ curves S. are constructed orthogonal 

to the C . The intersections of these curves with C define a set of 
n n 

points P. The boundary may now be approximated by the 

.traight line segments i, i+1 between the points P. ^ and n* 

coordinates (y. , z. ) of the points P. together with the corresponding 

n n * n * 

represent the basic input data which defines the surface geometry in 

the program. Denoting the number of segments in a cross section by iL 

and the number of stations by N the computations of this repo • have 

been carried out for N=10 and iL=20. 

From the points Pj jj intermediate points P/ ^ between 

P. and P., , on C are derived. It is assumed that the coordinates of 
n It i , n n 

P/ ^ may be represented by a Taylor's series in terms of the distance from 


^i, n» 


0U]Cl\-\j „ 
01' p: . 
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As = p 7 Tp 

i,n 


Reduction of this expression to one in terms of the discrete points ^ 
results in the following form (with a corresponding form for iz/ 

• . TV l(i. n) , ^^^i l(i, n)^ , 

where Dy. is obtained by first computing the divided difference 

^^i+i n " ^i taking this to represent dy/ds at the intermediate 

point p/ a distance l(i,n)/2 from P. . Linear interpolation of dy/ds 

between P/ ^ ^ and P/ ^ yields approximately dy/ds at P, ^ and this 

is denoted by Dy.. DDy. is the approximation to (d^y/ds^). determined 

by operating on Dy. in the same manner. In this report terms up to 

second order in l(i»n) have been employed. (Results obtained with P/ 

x^ n 

defined as above have been compared with those for Pf ^ defined simply 
as the mid point of the secant i,i+l and it was found that the latter case 
required double the number iL of segments to obtain comparable 
accuracy). 

21. Source density a 

a is determined by requiring that cp of Eq. (4) satisfy the boundary 
condition Eq. (2b) at a point P^ of each segment i,i+l of the boundary. 

The result of this process is a set of simultaneous equations for the den- 
sities O(i^n) at each segment i, i+1 of C^. There densities may be 
assigned to the pts. P/ » located as prescribed in Sect. 1. 

ip n 


iL 
i=i 


(^). = I A(j,i)a(i,n) 


(9) 
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where, referring to Fig. 3 for R(i,j,n) and 6(i, j,n) 

A0,i) * 2|8in[e(j.n) - 9(i, n)]ln[R(i+l, j. n)/R(i, j. n)] 

+ 6(i, j,n)cos[9(j,n) - 6{i,n)][ . 

The slope (dv/Sx)^ ^ may be written in terms of (^v^/Sx)^ ^ referred to 
body ^txes and the angles of attack 0. and sideslip ^ (see Fig. (4) ) 

(■^) ~ ® cos 6(j,n) - ¥ sin 9(j,n) . (10) 

j.n j»n 

Computation of (^v /^). from the surface data is described in 
o j, n 

Appendix A. 

Values of a(i, 10) obtained from Eq. (9) in the case of a circular 
cone at angle of attack are presented in Fig. (5). The analytical solution 
for 0 in the case of bodies of revolution is: 


= . 1 r s'(x)/2TTr , ^ 


2n 


cos 


This result is also presented in Fig. (5) for comparison. 

3. Potential cp 

Once the source density 0{i,n) is determined Eq. (4) yields an 
explicit representation of cp. Integrating over the segments i, i+1 of C^: 


iL 

co(j,n) = 2 ^ a(i, n)-JR-(i+l, j, n)* u(i, n)lnR(i+l, j, n) 
i=l 

- R(i, j,n)»u(i, n)lnR(i,j,n) - R(i, j, n)# n(i, n)6(i, j, n) + l(i, n 


= 2 


iL 

\ ■ 

La 

i=l 




( 11 ) 


Although w(j,n) is not of direct interest the auxiliary functions 

^(i. j. n)/CT(i, n) appear in the results for 5cp/3x and so must be computed. 
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4. Axial Potential Derivative Sq)/Sx 

dip/cfuc is obtained by differentiation' of Hue intdgta|iiafi Eq.E(^ to 
first obtain an exact expression which is then approximated by evaluating 
the result over die segmented boundary. This is felt to be preferable to 
the procedure of differentiating the approximation to cp given in Eq. (ll)t 
However, some care must be exercised when differentiating since the path 
of integration C(x) of the integral in Eq. (4) is itself a function of x. The 
details of this process are supplied in Appendix A. The resulting expres- 
sion for dcp/dx is found to be: 


- 2rJ ^ 


f (H) 

I A ° J 


in(Z-C)ds 


. X dC 

* 3 ^ 


which after integration over the segmented boundary yields: 


j, n j IL o Ji, n ' ’ ' 




The radius of curvature h(i^ n) and the derivatives (6a/6x)^, d^/ds, 

6v/5x are evaluated at the mid points of the segments i, i+1 by interpola- 
tion procedures described in Appendix A. 

Calculations of dcp/ ox for the circular cone at angle of attack are 
presented in Fig. (6). For comparison^ the analytical result for bodies 
of revolution is: 


9cp S*(x) , S^(x) _ o S(x) _ 

S “ TS“ -* — 


A plot of Eq. (14) for points cn the cone surface is also provided in Fig. (6). 
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S. Velocity Components v, w and q* = v* + 

}'" Differentiation of Eq. (4) with respect to Z yields the complex 
velocity function 


rc(c) 


ds 


(15) 


V - iw = - 2 0 

which, upon integration over the segmented boundary yields: 

v(j, n) - iw(j, n) = 2 o(i, n) [in + j, n)] (16) 

i 

is most conveniently found by noting that it is the sum of the squares of 
the normal and tangential velocity components. Thus, upon introducing the 
boundary condition Eq. (2b): 

q® = (■^) + (vcos6 + wsinQ)® (17) 

on the segment j, j+1 this becomes 


J.n 


2) a(i, n) 


i 


cos (e(j, n) - 9(i. n))ln 
- 6(i, j,n) sin(0(j, n) - 6(i,n))|| 


(18) 


6. Pressure Coefficient and g^(x) 

Cp depends upon q® and Sco/9x as determined above and the 
derivat g '(x). Differentiation of g(x) must be carried out with due 
concern for the nature of the improper integrals appearing in Eq. (5). The 
result of the differentiation process as given in Appendix A, Sect. (5) is; 

8'"n' = (^J " 

. suonnx^. S'(.))n(l-x„)} 


whe re 
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‘n* I I - s;,;ii.(jt^ - 


w ** — A 

Jn - J lif(V‘)®'<‘>* “ I (®;,M - s;;)l-<% - S^) 

c m=0 

*4, = l • 

m m+l m 

To compute the second derivatives of the cross sectional area required for 
g^(x) the first derivatives at are found by finite differences between 
Xjn a«d Second derivatives S'(x^) at x'(m) = 

are then found by finite differences between S* at x^ and Finally 

S^(x^) is determined by linear interpolation of between x^ and 


m+r 


Because of the possible singularity at x=0 the results are sensitive 


to the value of S*(0). Rather than compute this second derivative from 
discrete data it is assumed that the nose of the body may be specified 
analytically and that an analytically derived value is available for S'(0). 


The pressure coefficient 


c -^^=- 2 (^ + g'w)-q’ (20) 

p pU®/2 * 

may now be computed. The computational precision may be evaluated by 
comparison with the analytical results for a conical body of revoluttana- 
In this special case we obtain for points on the surface of the body 


g'(x) . . In (2,/*(I.x)) ^ 




+ 2a cos 9 + a® 


with 5cp/5x as given by Eq. (14). 
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Comp* 2 ted values of for the conical body are presented in Fig. (7) 
together with the analytical results obtained from Eq. (20) and Eqs. (14, 
21 . 22 ). 

7. Force and Moment Coefficients 

From Eqs. (7.8) for the force and moment coefficients it is seen 
that a determination of the "doublet" strength A^(x) is required. This 
term represents the coefficient of 1/Z in the expansion of the complex 
potential W(Z) about the origin (see Eq. (3) ). A^(x) as derived in 

Appendix A.Sect. (3) is given by: 




(23) 


where 




) . 


To obtain force and moment coefficients Aj(x) is substituted into Eqs. (7) 
and (8) which may now be written in a more convenient foirm by introducing' 


the centroid location Z in terms of body axes coordinates. 

g 

Z = Z - (ia + Y)x . 
g go 

The resulting force and moment equations are: 


(24) 


F + i F 

^ 5 = 2tt A,.(x) - (Y + ia)S + (Z S)' 

PU» 10' ' ' 'go 


(25) 


M +iM r F +iF ^ . 

—X. i-jx—i - I [2n Ajq(x) - (7 + ia)S]dx - Z . (26) 


PU® 


PU* 


Numerical evaluation of the integral in the expression for the moment 
coefficient is carried out by the trapezoidal rule using values of a^q(x^) 
and S(x^) obtained at each of the stations x^. Computation of Z^^(x)S(x) 
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is described in Appendix A* The derivative of Z S at x is obtained 

po n 

by first computing the divided difference between station? x^, and 
then letting this represent (^^S)^at x^. The derivative at is deter- 
mined by linear interpolation of between and . 

Analytical results in the case of bodies of revolution at angle of 
attack Ot are particularly simple: 


F 


— = aS(x) 

(27) 

PU® 


( \ 


— ^ = - a(x S(x)-V(x)j 

l28) 

pu^ ^ ' 



where V(x) is the volume of the body up to the station x. 

Computational results for the cone and ogive bodies of revolution 
at 01= • 1 are presented in Figs. 8 and 9, together with plots of Eqs. (27) 
and (28) for comparison. There results are presented in terms of the coef- 
ficients C- (x), and C, :(x) defined at the end of Section A. 

Lf M 
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C. Application to Typical Fttaelag[e 

A typical "si h-sided" fuselage toget)»jr ‘with>dQtaiIeiiaf><aiis$ of 
the geometry, is shown in Fig. 10. Cross-sections have been made of 
straight lines and circular arcs while the profile is composed of straight 
lines and parabolic arcs. Stations have been taken closer together ■ 
toward the rear of the body to p: omote a more accurate determination of 
total force and moment. Stations are situated farther apart over the center 
section since there is no change in cross -section for 1/3 < x <" 2/3. 

Processing of the surface data in accordance with paragraph 1 of 
section B is shown in Fig. 11. 

Results of the computation of pressure coe.ficient, force coefficient 
and moment coefficient are given in Figs. 12 and 13. 
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APPENDIX A 


DERIVATIONS 

1. Source Strength a 

Computation of 0(i, n) over the segment i,i+l proceeds by applying 
the boundary condition Eq. (2b) at each segment of C^. If ^cp = q = Jv + kw 
represents the velocity vector, the corresponding complex velocity in the 
crossflow plane is obtained by differentiation of W in Eq. (4) with respect 
to Z; 

C^(C)ds 


iw = - 2 




The contribution by the sources located on segment i,i+l to the velocity 

at P/ is first evaluated. Noting that i, i+1 makes an angle 0(i, n) with 

respect to the horizontal axis, we have 

jr j i®(i» n) 

dC = ds e ' ’ 

and the contribution to the integral in Eq. (Al) may be written: 


c. 


i+1, n 


^[v(j 


, n) -iw(j,n)]= - 2a(i, 


J 

^i, n 


, n 




(A2) 


After integration of the last term and summation over all contributing 
segments, the result may be written: 

v(j,n) - iw(j,n) = 2^a(i, n) ”^[ln (A3) 

i 

in which, referring to Fig. 3, the quantities R(i,j,n) and 6(i,j,n) are 
defined by the relationships: 

6{i,j,n) = i '(i,j,n) - 'I'(i,j, n) . 

To insure uniqueness of the complex velocity, care must be 
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exercised in assigning values to the angles ir(i, j, n) and t|r'(i, j,n). Refer- 
ring to Fig. 3, these are measured counter-clockwise from the positive 

y axis so that when facing from P. to P.,, , a point P/ just to the 

^ ® i,n x+i,n ^ 

left of i, i+1 shall define an angle ^'(i, j,n) = 9(i,n). As p/ traverses 

J»n 

a path around P. to a point just to the right of i,i+l, iKi» increases 
1 , n 

from 0(i, n) to 9(i, n) + Zn. The same holds true for i!r^(i,j,n) as P/ 
traverses a path aroung consequence of these definitions 6(i,j,n) 

becomes -tt when approaching i,i+l from the right and rr when approaching 
from the left. This discontinuity reflects that exhibited by the stream func- 
tion upon traversing any closed path which encloses a distribution of finite 
sources. 

From the boundary condition Eq. (2b), we have: 

(•^) = v(j, n) sin 0(j, n) - w(j, n) cos 0(j, n) . (A4) 

j.n 

After substitution of v and w from Eq. (A3), this last expression becomes 
(|^) =^a(j,i)a(i.n) (A5) 

j,n i 


where 

a(j,i) = 2{sin(0(j,n) - 0(i, n))ln 

+ 6(i,j,n)cos(0(j,n) - 0(i, n)j| 

In addition, we see from Fig. 4 that the slope 3v /2x may be expresced in 
terms of the body slope 5v^/9x referred to body axes: 


(■^) = ( 15 ^) + acos 0(j, n) - Y sin0(j,n) (A6) 

j, n j.n 

thus eliminating the necessity of constructing a new set of projections 
similar to Fig. 2 for each set of a and Y. Satisfying Eq. (A5) at each of the 
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points P/ on a given cross-sectional boundary yields a set of equations 
for a(i,n). 

2. Determinatijn of q», 9q)/dx 

A knowledge of o{l, n) allows the numerical integration of Eq. (4) 
for o in a manner similar to that for the complex velocity above: 


^i+l,n 

®(j. n» = - 2Re^ c(i. n) j n “ 

i C. _ 


(A7) 


After integration, cn(j,n) may be written concisely in the nomenclature of 
Fig. 3: 

V f- 

?!(j, n) = 2^ o(i, n) |R(i+l, j. n)* u(i, n) In R(i+1, j, n) 


- R(i, j, n)- u(i, n)lnR(i, j, n) - R(i, j, n)* n(i, n) 6 (i, j, n) 


f l(i, n)} = 2 y 0(i. n) (A8) 

i 

in which use has been made of the geometric relationship: 

R(i. j. «)* n(i, n) = ^(i+1 , j, n)* n(i, n) . 

The derivation of &on/9x must take into account the fact that the 
path of integration in Eq, (4) is a function of x. Referring to Fig. 1, we 
shall distinguish between increments of a dependent variable taken along 
C(x) and denoted by d( ) and increm-^nts taken normal to C and denoted by 
6( ), Differentiation of Eq. (4) then yields 


5 ^! 


- 2Re{<f |^ln(Z - C)ds - <f ^ 


0(C) K 

C Sx 


ds 


+ ^ 0(C) Ir (Z - C) 

From Fig. 1 it becomes evident that 


(A9) 
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f fd»> • 0 v'tiO = 6,>i ^ (A10> 

^hei« h(C) i« tk» t^itts of curvatur* of C(x) at Q. In addition, w hare 
from Fig. I, 

|1 , ,i (9 - n/2) 

To tvaluate Ccf/6x we motie, reforring to Fife. 4, 

(A12) 


64 - 4_-0(i,n) 

37 - 6k * 

6x-H) 


wheae 0" de notea the value of at the point P*'. The relative displaca- 

maat between F. _ and P* is shown in Fig. 4, as !t would appear in "wind" 
1, n 

* 

axoa. However, the cempukatioa of d has been carried out in a body axis 
frame af rsieyence. To make use of the results of that cotuputation we note 

it 9 

•lat O im the wind axis frame corresponds to d in the body ax*j frame. 
From Fig. 4 than* we have 


= d(i,n+l) + ^ (01 sin 9 + Y cos i)6x 


ds 


(A 1 3) 


which, aftea substitution into Eq. (A 12)^ leads to the required expression, 

(A14) 


If = (te tf^iastaetfeose! 


‘67' ■ ds 

o 


where (-^) is the dertvative evaluated in the body axis frame. Finally, 
o 

intaoducing Eqs. (AlO), (All), and (A14) into Eq. (A9), 

+ff(*8in8 + ¥cos8 + j|^]ln(Z-C)de 

n 

Again, assurnlng that quantities in the brackets of the integrands are 
Fonstant over i*!4 1, the integrations proceed in a straightforward manner: 


origin AT. T’A' -. H ^ 

OF FOOR UUiU-.U* 
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(AU.) 

i , n I 

. 4 

in which we note that (6v/6x) s ?v/?.\ as defined in Eq. (A6). 

Equations defining (dc/ds), (6.7/?x) , fv /6x and 1/h at the 

o o 

point P/ ^ are provided in Sections C-1 and C-3 in Part II of this report. 

A description of computational process is given here: 

a) cro/ds - CT at P. is first obtained by interpolation be twee .1 

1 , n ' r 

the computed values of cr(i^ n) at da/ds at P/ ^ is then set equal 

to the divided difference between these ^'nterpolated values of CT. (see 
Section cfis of Part II). 

b) (C(?/6 k)^ “ the derivative at the mid-poir. of the interval 

*n* ^ntl equal to the divided difference between rr(i, n) and c(i^n+l). 

Linear interpolation between these derivatives then yields (5c/^x)^ at x^. 
(see Fig. 14 and Section C-3 of Part II). 


c) 6v^/fx - Referring to Fig. 15, the displacement is deter- 
mined by interpolation bet^'een ^ C. and ^ T.,. . 6r/(x - x ) then 

^ ^ w iTl,n nti n 

represents /6x at Interpolaftoti between the stations then 

yields fv /f»x at x (see Section C-1 of Part II). 
on 

d) T/h - 6 at ^ is determined by interpolation between values 

of ®(t, n) at P/ . The curvature 1/h at P/ ts then set equal to the 
* 1 , n 1 , n 

divided difference between f at P. , and 6 at P. . (see Lection C-3 

11 ^ n 1 , n 


of Part 11). 


<•) 
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3. "Doublet Strength" A^(x) 

Aj(x) is the coefficient of the 1/Z term in theexpansion of the com- 
plex potential W(Z1 about the origin (see Eq. (3) ). If the integral repre- 
sentation of W from Eq. (4) is expanded we find: 

W(Z) = (-2| ff(C)ds)lnZ + (Z| C tr(C)ds)i 

+ (ZX C* c(C)ds)-L + (A17) 

z® 

Thu8p we have for the coefficient of the 1/Z term: 

Aj(x) = 2 1 C o{Q)ds 
Introducing body axes coordinates 

C • (ia + T)x 


we have 

Aj(x) = 2(p cr(C^)ds - 2(ia + Y) <y(C^)ds 


The last integral on the right hand side is recognized as the coefficient of 
the ’‘source” term in the above expansion of W(Z). According to slender 
body theory Ref.(l), this is related to the rate of change of cross-sectional 
area: 


2 ^ a{Q^)da 


S^(x) 

2rr 


our final expression for the ^doublet” term is therefore 

Aj(x) = 2^ Cq + (ia + P)x 

Integrating over the segmented boundary C^. 

i+1 

f cr(Co)ds = Y a(i, n) J Cq ds 
i 1 

the last integral may be interpreted as the mome'- ' itl.e . rc i, i+1 about 
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(AI8) 


the origin and may be approximated by (y/ -f iz! ) l(i, n) so that 

I, n n 

Aj(x) = 2^^ a(i,n) l(i,n) (y^ „) 

4. Cross-sectional Properties 

Computation of S(x^), Zg^S(x^) and their derivatives is accomplished 


with the aid of Stokes' theorem in the complex plane. Thus, 

JT T ‘o°^o 

** “ (A20) 


S(x) = dC^ 


2 S(x) = o C-T^dC 

so j o o ( 


go ' ■ 21 j ”o “o "o 

which expressions, after integration around C^, yield 

S{x ) = 4 y (y/ dz. _ - a/ „ dy. ) 

n 2^ 'i,n i,n i,n i,n 


and 


1^ 




go ^ ^ 1 , n 1 , n 


(A21) 


(A22) 


where 

r®. = (y/ f + (z/ )« 

i,n 'i, n 1 , n 

“ 2 « 

i,n 1 + 1 , n 1 , n 

*^^1, n ” ^i+I, n ” ^i, n 

5. g ^(x) 

The derivative of g(x) appears in the expression for the local pres, 
surf coefficient, Eq. (6). To avoid the occurrence of singular integrals, 
differentiation is accomplished by first integrating by parts the integrals 
appearing in Eq. (5) for g(x) and then differentiating the resulting expres- 
sions 

X X 

f S'(t) In (x-t)dt = - S"(0) (x-xlnx) - fs'"(t)[(x-t) - (x-t) In (x-t)]dt 


'oUALVl-'^ 
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then 


S'(t)ln(x-t)dt s S'(0)lnx + f S"(t) In (x-t)dt 


and similarly 
1 


r S'(t)ln(t-x)dts - S'(l)ln(Ux) + [ S*'(t) In (t-x)dt 

X X 

Thus, differentiation of Eq. (5) for g(x) yields: 

1 

g '(x) = ^ {s^x) In (|) + I f S*(t) In (t-x)dt 


- I J S"(t) In (x-t)dt - * i 

o 

. s'(i) 1 s'(0K s'(iK A 

Expressing the integrals as Stieltjes integrals facilitates their computation. 
1 N-1 

‘n = f S’(t) = Y . S^yin te' - *_,) 


and 


n 


m=n 


n-1 


= j l„(x„.t)dS'(t) = £ - s;;)ln(x„ - x^) 


m=0 


where x' = (x +x ,,)/2 
m m m+1 

we thus have 


n n 

-S"(0)lnx - S'{l)ln(l.x )\ 

n n j 


(A26) 


The occurence of singularities in g(x) and g '(x) at x=0, 1 signifies the 
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failure of slender body theory in these regions unless S is sufficiently 
well behaved there i.e-., first and second derivatives equal to 0. For pointed 
bodies S'(0) = 0 and the occurrence of S'(l) s 0 is common. 
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PART n 

FORTRAN PROGRAM 

A. Input 

1. Comments 

The body axes coordinates y. , z. at x may be read from 

n* X, n n ' 

cards or computed by a code supplied by the user; the indices IX and 

IR are set equal to 0 or 1 depending upon the choice made. After the 

source strength c is computed the program computes CD, 3co/3x, at 

the locations P/ ^ on the surface. The capability of computing these 

qtxantities at arbitrarily specified points on or off the body has also been 

included to facilitate induced flow studies. Thus cn, 3co/3x, C are 

P 

computed at P/ ^ or at locations supplied by the user as additional input, 
depending upon whether the index lYPP is set equal to 0 or 1 

2. List of Fortran Symbols for Input Data 


ALP 

BET 

ACH 

SPPO 


SREF 

ENG 

REFL 


Angle of attack a, positive for nose up attitude relative 
to wind axes. 

Angle of yaw positive for clockwise rotation 
about z>axis. 

Free Stream Mach No. 

S'(0) Second derivative of cross-section area evaluated 
at the nose. It is assumed that this is available from 
analytical considerations regarding the special geo- 
metry of the nose section. 

Dimensional reference area. 

Dimensional body length. 

Dimensional reference length. 
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lYPP 

IL 

NL 

IR 


IX 


>0 if coordinates of p/ ^ are compoted by program, 

=1 if P*. are to be read from ioout cards. 

1, a ^ 

Number of segments into which a cross-sectional 
boundary is divided . 

Number of longitudinal stations at which cross-sections 
are taken. 

=1 if y- . z. are to be read from input cards. 

=0 if these cards are to be computed hy a code 
inserted after statement 111. 

=1 if are to be read from input cards. 

=0 if these stations are to be computed by a codd'te' 
inserted after statement 113. 


ISYMLR = 0 
= 1 

ISYMUD = 0 
= 1 


if contour does not have lateral symmetry 
if contour has lateral symmetry 
if contour does not have vertical symmetry 
if contour has vertical symmetry 


ISR 

X(N) 

Y(I.N) 

Z(I,N) 

YPP(D 

ZPP(I) 

3. 


if ^ 0 SREF will be defined = S(ISR) 

Dimensional longitudinal coordinates x^. 

Dimensional coordinate y. 

i,n 

Dimensional coordinate z. 

i,n 

Dimensional coordinate of collocation pt. 
Dimensional coordinate of collocation pt. 
Preparation of Input Cards 


/ 

z. 

i,n 


Card # 


Format 


Variable 


1 5E15.8 


2 5E15.8 


ALP 

BET 

ACH 

SPPO 

SREF 

ENG 

REFL 
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3 


1015 


lYPP 
IL 
NL 

m 

IX 

ISYMLAR 
ISYMUD 
ISR 

The following cards are prepared in the order presented, when the indices 

DC, IR, lYPP are as specified 

If IX=1 10F8.0 X(l) 

X(2) 


X*(NL) 


If IR=1 10F8.0 

If ISMLR = 1, ISYMUD = 0, or 1 
1= 1 placed in 4th quadrant 
I = IL placed in 3rd quadrant 

If ISYMLR = 0 ISYMUD = 1 

I s 1 placed in 1st quadrant 
I s IL placed in 4th quadrant 


Y(l,l) 

Y(l,2) 


Y(I.NL) 
2 ( 1 , 1 ) 
2 ( 1 , 2 ) 
2(1, NL) 
Y(2, 1) 
'^( 2 , 1 ) 


If ISYMLR = ISYMUD s 0 

no restriction on placement of 1=1 Y(2, NL) 

2 ( 2 , 1 ) 


2(2, NL) 


2*(IL,NL) 

If lYPP = 1 5E15. 8 YPP(l) 

2PP(1) 

YPP(2) 

2PP(2) 


y’pp(il) 

2» !(IL) 


If IR=0 A code to compute ^ must be inserted after 

statement 111. 

if 1X=0 i'. code to compute must be inserted after statement 

113. 
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B. Otttpttt 

1« Input parameters 

The first row of output presents the pertinent input parameters 
ALPHA. BETA, MACH NO. , SPP(O), REF AREA, BODY LENGTH, REF 
LENGTH. 

0(j,n) and Cfi(j,n) at the location y/ 2 / are presented as 

J»n J.“ ^ 

follows for 1 < n < N 


SIGMA 

a(l,n) 

a(8,n) - - -ffOL.n) 

PHI 

cna,n) 

orj(8,n) - - - cn(IL,n) 

Y PRIME 
Z PRIME 




* - - ^IL.n 
■ " “ *'lL,n 


0(7, n) 


®(7,n) 



oinr.TN-M. 
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3 . 5eo/2bc 


(9cp/5x). 

'j, n 


at the points 



are presented as follows: 


D Pm/D X 


(dcp/dx)^ J 

(S!p/Sx)g’, .... , 


(5q>/Sx)j 


4. ARjq(X^) , AIjq(X^) 

Rea] and imaginary parts of the "doublet strength" Ajq(x^) are 
presented as follows: 

AI AND AR 

A1jq(Xj), AR^q(xj), ^jq(x2)» ARjq(x2) - - - AIjq(x^) 

ARiq(x4) - -- -- -- - 

5. Force and Moment coefficients, treasure Coefficient 

Pressure coefficient C at P.^ is computed for 1 < n < N-1. 

p j,n ^ - - 

Force and moment coefficients are presented as follows: 

N = n, CY = Cy(x^), CL = Cj^(x^), CN = C^(x-) 

Vn^» GP=gX) 


Cp(l,n) Cp(7,n) 

C (8,n) ----- C (IL,n) 

P P 
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C. Sommary of Programmed Equations 

Thaaa equations are presented in order of use. The Fortran tvmbol 
at the left represents the quantity at the left hand side of each equation. 

1) Computation of o(i, n) 


^ILP.N) 

yiL+i,n*yi,n 


Z(1LP, N) 

*iL+l,n s s. „ 
i,n 


Y(ILP. N) 

nL+1 * ^2.0 


Y(IL2, N) 

yiL+2 = y2.n 


Y0L3.N) 

^iL+3 “ ^3, n 


Y(IL4. N) 

^iL+3 * ^4. n 


FKILP, N) 

l(iL+l) = l(l,n) 


FK1L2, N) 

KiL+2) s 1(2. n) 


FKIL3, N) 

l(iL+3) = 1(3, n) 


DPY(I) 

“'n = <n+l.n - 

} < i < iL ^ 3 

DY(I) 

D Vj. iKi» n) + D Vi * ») 

Dy. * ^ — 

' l(i,n)+ l(i-l,n) 

2 <i 3 

DPY(I) 

DVi* (Dyi+i - Dy.)/l(i,n) 

2 < i.S 2 

YP(I) 

DVi jl(i.n) +DVil(i-l.“) 
l(i,n) +1(T-T7nl 

3 < i < iL + 2 

YP(1) 


. 3 < i < iL + 2 

yp<i) 

I'l ' Vl+i 


YP(2) 

’'i * n'L+2 
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The abowe operittions from Y(ILP, N) to YP(2) are repeated for Z(ILP, N) 
to ZP(2) to obtain z/. 


R(I,J) 

R(i,j.n) =f(y/ -y. + (z/ - 

•’ I'j.n 'i.n' j,n 

FL(I,N) 

*«-"> ■ ’'i.n* * ■ 

ST (I) 


CT(I) 

co»9(i,n) = (yj^ - y )/I(i,n) 


For the computation of angles it is assumed that a computer Will obey the 
following rules: 


0 

< sin*^ sin 6 < .t/2 , 

sin 9 (4) 

-tt/2 

< sin* ^ sin 6 < 0 , 

sin 9 (-) 

0 

<cos*^cos9 <rr/2 , 

cos 9 (+) 

tt/2 

< cos* ^ cos 0 < n , 

cos 9 (-) 


T(I, N) 


AS 

ZZ 

YY 


9(i, n) =< 


sin”^ sin9(i, n) 
TT-sin”^ sin9(i, n) 
-sin” ^ sin 6(i, n) 
2n-sin” ^ sin 9(i, n) 


sinv(i,j,n) = (z/ - z )/R(i,j,n) 

j, n ^ 

Az = z/ - z. 

j,n i,n 

Ay = y/ - y. 

J.n i.n 


sin 9(i, n) 

cos 9(», n) 

> 0 

> 

— 

■”* 

> 0 

<0 

< 0 

< 0 

<0 

>0 



Az 

Ay 

sin*^ sinY(i, j, n) 

> 0 

> 0 

TT-sin* ^ sin y(i, j, n) 

> 0 

< 0 

. .1 . . 
TT-sxn 8inY(i,j,n 

< 0 

< 0 

2TT+sin* ^ sin Y(i, j, n) 

<0 1 

> 0 
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♦ (ij.n)* 


PHS 


t'(i,j.n) = 


Y(i,j.n) 
v(i, j.n) + 2 

v(i+l J,n) 

*(i» j»n) 


Y(i,j.n) > 6(i, n) 
v(i,j,n) < 0(i, n) 


Y(i+l,j,n) > 0(i,n) 
v(i+l J,n) = 0(i, n) 


v(i+l,j,n) + 2 tt , v(i+l,j, n) < 9(i, n) 


D(J,I,N) 6(iJ,n)= n) - ❖(ij.n) , i ^ j 

D(J,I,N) e(j.j.n)s-TT 

The following redefinitions of 0(i, n) assure continuity of 0(i, n) when 
passing directly between first and fourth quadrants: 

0(iL+l,n)s 0(1, n) 

A0 = 0(i+l,n) - 9{i,n) 

0(i+l, n) + 2n , A0 < -(f + 10"®) 

0(i+l,n)= ^0(i+I,n) , - f < &6 < f 

0(i+l,n) - 2 f , A0 > F + lO""* 


FL(IL+1, N) l(iL+l, n) = i(i, n) 


BE(I,N) 

6E(1. N) 
DR(I) 


DNX(I) 


K(itl . n) = 8 ( 1 . n) f 


h( 1, n) - >t(iL-fl,n) - 2^ 


6''„li.n) = (Vj sinx(i,n) 



1 < n < n-1 

fCv ] + 5v^(i+l,n)]/2 

\3x/ “ X . , - X 

n+1 n 

2 <n < N-l 
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IMII.N) 




2 < n < N 1 


, Sv 


i*W i,N-l 


PN(1,1> 



(I*) "(is?) ♦ n) - ^ sinfi(i,n) 

i.n 

3in(ff<j,n) - #^i,n)l = si«#(j,n) cor 6(i, n) - sin6(i,n) cos6(j,n) 
co«[9(), n) - n)| = cosl(j,n) cos ^’,n) + sin6(j,n) sin6(i,n) 

a{i,j,n) 2 ■[•4f«(8(j, n) - »(t, n)] In 

+ co8[f<j, n) - 6(i, »)] 6{i, j, n) 

0(i,n,|= 1 

,i 

2) Cotti|^tation of cnfi^n) 

RT lR(i, j,n)*u(i,n| = (y/ - Yj cosS(i,n)+ „) sin«(i,n) 

, R« R(i^^ j,n>*u(i,n) = (y.' ^ - y.^ ^)cos6(i,n) 

■*■ « * *i*i J *»«•(*. n) 

j, n It 1 4 

RN Rd*, j,n).n(i, n> = - (yf _ - Y. J n) 

.144^ ^ 

* n ■ J 

3 , n 4* « 

DTdJ.I.N) {-^^1 * R(i+l,j,n)*u(i,n) InR^iYl, j,«) 

‘ R(i,j,n)*u{i,n) lnR{I,j,n) 

• j, n)* n{i, n) 5 (i, i, ») ♦ 4{-i, «) 


DN(t,«) 

rrT 

CTT 

Am ) 

StCdl.Nf 


u 


# 


mm mm - 



PH(J) 


ORIGINAL' PAGE IS 
OF POOR QUALITY 




3) Computation of (do)/3x). 


SIGOLP, N) a(iL + 1, n) = o(l, n) 


ss 


DSX 


DoX 


RD 


TX(IN) 


fdo\ _ '^(i+l.n) - ^i.n) 


XIN(I) (^) = g(^» "tlL^ ") 

' 413C f • XI - m • ;X 


n+l 


r^o) . /Ao\ (La\ 1 n ‘n-l 

\^x/ . ” '"Sx/. . l^Axi, " \T^A |Jx . - X I 

oi, n i, n-1 i, n i, n-1 n+l n-1 


(If). =(^). . 

oi, n i,n-l 

1 tuiz - Mi+l,n) - H(i.n) 
l/h(i,n) 

iLl 


i.n 


S). = 4 (asine(i,n) +f cose(i,n))(g) 

■ “ i i o». « li 

a(i,n)^6v'i »/'6v', e,. . 


i.n 


1. n 


4) Computation of q^(j.n) 
Q2( J, N) 


(j»n) = +I 2 ,'' a(i, nlFcosfefj.n) - 6(i,n)) 

j,n ( i '• 

«in(0(j,n) - 6(i, n)jj 


2<n< N-I 


2<n<N-l 


n) 
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5) Computation of Crosa -sectional Properties 


Y2P 


a (y/ )* + (a/ )* 

''i, n' i, n' 

ZZ 

dz. 

x,n 

* *i+l, n ” *i, n 

YY 


^i+1, n “ n 

S(N) 

S(x„) 

(y.* d*. - z! dy. )fZ 

L. 'i,n J,n j, n 'i, n' 

i 

SYG(N) 

y S(x ) 
'go n' 

dz. 

. i.n 1 , n 

SZA(N) 

go n 

- * V r» 

dy. 

. \,n 'i.n 

DSYG 

DYS 

n 

s JK2 StL: 82 !!_ , l<n<N-l 

*n+l ■ '‘n 

SYP 


= DYS„_, ♦fDYS„-DYS„.,)^“ // 

n+1 n-1 



2 < n < N-1 

SYP 

(v<*. 0 ' 

« / 

a 2 DYSj - (yj^S(2) j 

SYP 


= 2 DYSj,., - (ygS(N-l))' 


repeat for 


SPXP 

XP(J) 

SPPXPP(J) 


S'(x' ) 
m 


# 



S'(x') 

m 



- S(x ) 
' m 

^m+1 

- X ^ 

m 

<’‘m+l * 


S'(xj') - 

S(Xj)/xj 

/ 

’‘l - 

Xj/2 

S'(x' ) . 
ni 

• S'(x' ,) 

' m- 1 


X - X - 

m m- 1 


1 <m < N-1 
1 < m < N- 1 


SPPXPP 


2 < m < N- 1 



XPP(J) 


X* =s (x * I + » 

m m-1 m * 


OK‘IGIXAL PAGE IS 
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XPP(l) 

SPPX(J) 


X, = 


S* 

m 


(Xj + x^/2)/2 


m ' « » 

X . , - X 

m +1 m 


1 < m < N-2 


SPPX(NM) 


SPPX(NL) 


SPX 


*N-1 






*N-1 ■ * N-2 


■ ^N-1 ^ f^N-1 ■ ^N-2^ 


^N-1^ 
*N-1 ' *N-2 


s'(i) = s^(x^ ^) [s^(x^ - s^(xj;^ ^)) 


’^N-l ■ ^^N-2 


6) Computation of g^(x), C , ^*{ 0 ) assumed = 0) 

P — 

N-1 


RIN 


n m + 1 m m n 


m=n 


n-l 


RJN 


GP 


n ^ ititI m n m 
m=0 


n 

-S^(0)lnx . S^(l)ln(l-x )\ 

n n ’ 


1 < n < N-l 


Wl(J) 


C (j.n) = - 2 (^) - q»0,n) - 2gV^) 

P 1 , n 
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7) Computation of Force and Moment Coefficients 
AR(N) ARjq(x^) s 2 ^ o(i, n) 1 (i, n) ^ 


AI(N) 




W3 

we 

W2 

W3 


Fy/pU» = 2 ttARjq(x^) - T S(x J + (y_ S(x J ' 


n -"go ' n 


F,/PU» 


( N 9 

% S(x ) ) 
go n'/ 




C = 2(F /pU=)(L»/S^^,) 


SUM f[2nARjQ(x) -Y S(x)]dx = 2 ttARjq(Xj) Xj/2 


?'V 

X Zn I }[aRio(x^^j) - Y S(x^^j)/2rj 
m=l 


X 

SUMl I* [2tt AIjq(x) - a S(x)]dx=[2TTAIjQ(Xj) Xj/2 _XjaS(Xj)/2] 

o 

m=r 




•’‘mtl - 
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W5 


W4 


W5 

W4 


page is 

X ' IVuit QUALITY 

My/pU»= - x(F^/pU*) + J[2 ttAIjq(x) - a S(x)]dx + S(x^) 


X 

M^/pU®= x(Fy/pU* - r [ZttARjq(x) - P S(x)]dx - S(x^) 

o 


Cm = S^,) 

C„=-2(M^/pU»)(lVL^jS„,) 
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D* Program Listing 



*• 

- 



0001 

’ 


oiur^sioN r (.i(i«i6 i tZ( '*c« 10 ) (ONfiiuiM »/Lrio*io) tZf'OO)* 

00006000 



tnS(30*JO) »AJ1 (^00) *SI C3U) t 

ooooooio 



(30) (30o jO) «r>(30«j )• it>) ,,»r IJ(;* J0*1M f sr (30) •X(l 6 ) 9 

ooocoo^o 




Wh( 1 -)) .T 1 n)*lr') tGS»iM( ^u*30) *Cr (30) • 

oooauojo 



4TX { ( js 1^) slfK (30) *44 (Ifi) 9ftr 1 16) *0.J<( 30) ••JE ( JO9 Inl tSC i^) • 

oonooo4(j 



^Sv( ( (^^1 /C»( 10 309 161 

oaooou^u 

000? 



A9*( IM • do) 9V*>^»A*>P do) .SPPA ( 16) 

00f)fl00n0 

000.1 



CI^tVSlON C(3»>) *M( )0) *‘)W‘>(30) •41(30)9X16(30) 

OOOOOU70 

0004 



Olcr:6STt> IW (30) tOV (3in •OP2(30I «wZ ('^0) 

OOnOOOt^O 

OOOS 



r*Ifc*f«v<;iOvi O.jTCl) 

OOHOdOnO 

ooo^ 



Ewbl vA(,r <CE (AjI.OUT) 

oooooloo 

0007 



EdUlV^Lf 6Ce (P«kSI • (P9PS) • (0«^*0( 1*1 tl6) ) • (GSUNtOT ( 1 • 1 • l6) ) • 

oooonuo 




1 (Ov tCT) 9 (O^aSi) 9 (.)»^*SI *^*bla( 1 • 14) ) 9 (^H9 - 1 *L) * (*»hY .i,SUM ( 1 9 1 ) ) * 

00(i00li?0 



. -»• (l« 1) ) • (1 *^) ) « < U 3) ) • < 1*4) ) * 

00000130 



3(0tw*>.'V/>( J • lo) ) 9 *( l*#^) ) 9 (ONX^AlN) 

OOOOOUU 


C 

oo» 

If trt-P^Q 4 ill i^t SET = To VP ZPPa2P hsPS 0=OS a60 P 

rPOOOOOlSa 


c 

««« 

[F lVPP=i T?tE6 VPP A^n 1^0 Mu'iT ►HE input 

00000160 


c 


IF 1^ = 0 T»,ln ^'OO-: to ro*-<PuTE r*Z “*iST «E iNSEPTfO aetfr 

ooooor/0 


c. 

«O40 

STaTP .FNT ill 

onooo XbO 


c 

«»«■ 

IF 1-^1. Tr*t^^ y*Z ''UsT Input 

OOOOf'IgO 


C 


IF T^tN COOF ro ro^POTE X must me INSERTlO aftew 

OOOOuPuO 


C 


113 

00000^10 


c 

«09 

I^ Ia- 1 TrttN X Mfi^t mE T mmuT 

00000^*4^0 


c 


IE IST'^u^-O TmI-ME n NO L*-FT TO wIGHl i>YMHETRY 

OOOOOP3J 


c 

ft 

IF (^>' •{>31 T^FME is left to ->15^1 sVmmET^^V 

oooaap-9u 


c 

ftftft 

IF I = TmE»*E l4 NO i>p TO 'JOft * svf 3£.fgy 

onoooMsu 


c 

ftftft 

IF Is>**OI»rl !<; ire TO 00*N >VMMPTr^V 

oooor 


c 


If • vLF- :n.'\JT VmLJF 

unofioM/y 


i: 

ftftft 

IE nOT = 'j ^-l«F >IlL -c RFCEETvFl) =S(IsR) 

oor.ooM^o 

000^4 



PE4.)('9^'*k*l) ff*iC*>9'»9'Ki;9S'^r, F*r‘N(>9 kEFL 

onooo<*-^i* 

OOOO 



*PI TE (4 9'>V'*) Ai T9if:H9SPM0 9SMf‘F9E hUwFFL 

nooGO iuo 

00l«> 


sog 

K(>Mva 7 r Vh ALMf^Ar *F 4 . 1 9^^ fA^9E>», W-jAs >*h».w*Lh NU*«*F6*<i9/ 

onroo310 




1 7hSh*^ (ni'tFir'-^r**-)^* •<•-*^6 E‘ AMEA = *F»^«)9*^^91 /'MnOOV Lt N(» 1 M = 9F6 *394X9 

00 non <<rii 




l-^EF L^•v^>^“'=•F^•39//) 

OOCOO ^^o 

OOil 



mE;. WS94f<0) IY*^k9 IL 9\L 9 In9 1X9 T SY' LU* I SY'^uf)# ISM 

oonoo ><»o 

oni? 



MI = T Au( 1 ,r.) 

UOOOO i'DiJ 

0013 




joroo <so 

00|<o 



e • I .r *3 

oooon ( /y 

00l«' 




00000 0 

OOM 



Nl '= iL-l 


001 1 



- I *• -•• • I 

f»0 0 0 0** II 0 

001- 



IF ( I >Y ’LH *t .(• ». ANO. r»YMM').FO. •!) ILL = IL 

oijnonH 1 0 

001^^ 



IE (I ) ILL = TL//^ 

GO', )0‘*c:0 

00^0 



IF f I ^ .E\J . 1 • «N 1. ( * >«r >• 'W ! L' •= 1 L/^ 

oOGi'U;'* 10 

0«VJ 



(c ( f ^ / i.Eij.n rLL-iL/4 

U00M0*9^0 

00/* 



IE ( r f *0) <ir TO 1 1 4 

OG-Yon^^^O 

oo^^ 



eE *!■)( ^ r « (n) tM-: 1 9NL) 

OGOO()-*r 0 

00^4 



«-r To n- 

OOOnoH/O 



1 1 1 

Cl)'" T I %r tf 

») a 


c 

ft ft ft 

IE 4 -.*»S :>f.n Tei .* COr)f TO CO^'^UIE MUST hE 



c 

ft ftft 

I )\.s^ -7 fk. 0 r*t iVE 

n 0 0 r 0 s j 0 

00/*^ 


I 14 

IF ( i-f .» VJ..1) 'h' ro 111 

onf'UGsi y 

0O?7 



no I 1 9N{, 

GOOOOS^'U 



lo^ 

r- , 0 ( 1 ( > f t • M ♦ I - 1 • I LL 1 

ooGons (y 




) ;0 1 1/^ ^ ' 1 • "ll. 

fj fl 0 0 t)S « 0 

0030 


1 .^6 

» eE ,1.' ( -S) ( / M .N) 9 r 5l 9 M L ) 

oiuroosso 
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. .'K QUALITY 


00:^1 

SOS 


OllOOOSoO 

003^ 


CtO to U? 

OOO00S7O 

0033 

m 

continue 

OOOOOSbO 


C '»44 

ir Y 4.NI) 7 «e>K NOT INPUT TrtKN CODE TO CON»^UTE MUST HE 

OOOOOS^iO 


C *44 

iNSE-fTKD HEHfc 

00000600 

0034 

U? 

wO b-^7 N=1 *NL 

ooooo^lo* 

003S 


^ (.i) (N) /tNG 

ooooo6;f0 

00 3n 


DO ^>^^7 I«l*TLL 

O00C063O 

0037 


V(l*N)=Yn*N)/ENO 

00000640 

0038 

SS7 

/(I.N)5/CItN)/KNG 

00000650 

003«i 


IFMLL.FQ^lL) GO TO 461 

OOOOObbO 

0040 


IFClLL.Nt.lL/4) GO TO 4SS 

00000670 

0041 


Irt=lLfl 

00000660 

004? 


ILL=TL/? 

000006«^0 

0043 


DO ‘*50 Nsl»NL 

00000700 

0044 


CO I = lb«lLL 

00000710 

004S 


i»>=n.L^i-i 

0000U720 

00«»#> 


V(f.M=VClO*N) 

00000730 

0047 

4S0 

/CI*>.)=-^IU>*N) 

00000740 

004M 

4SS 

Ir.= lLL*l 

000007b0 

004«3 


FYs-1 . 

O00UO760 

ooso 


F?=l. 

00000770 

OObl 


IFUSyMLH.KG.I) go TO 46b 

00000760 

005? 


FV = l. 

o»>noo7vo 

00S3 



ooooosuo 

00S4 

4bS 

DO N = )*NL 

OOOOOMIO 

ooss 


00 :*bO I = lJ^*Tl 

00000«4i0 

oos^ 



COOOObJO 

00S7 


V ( I .N) =Y ( Iu*N|4F y 

0OOOOH-*O 

OOS^ 

4&0 


00O006SU 

oos^^ 

4f^l 

roNTiNt't 

00000660 

0000 

1 

no US NsUNL 

00000670 

OOM 


Ak C ^r) =0.0 

oooooaMO 


\ 

AI ( ') =0.0 

OOOOObVO 

OObl 

\ 

S ( . ) = 0 . 0 

00000^00 

0004. 

\ 

S/ =0.0 

OOOOONlO 

00n«- 


SV.,(\)=').0 

OOOOO'^cO 

OOo^ 


IF C iTM.J.FO.l ) ht A>MS*b.>n 1 ) UPHU ) « t = l • IL) 

00000^30 

0047 


Y ( =Y ( 1 .N) 

0f)000'#«40 



/ ( lL^*v*)=/{i*N) 

00000*^50 

004«f 


!»• 1) GO ro i 

OOOOO'^bO 

00 70 


=< 

onnoO'^TO 

oon 


Nr =l 

oni»oo'>rto 

oor? 


*0 TO b 

OOOOO^NO 

0073 

3 

IF ( i.Ne .NL) t’O TO 

00*10) 000 

00 74 


Nl = >* 

ooooiolo 

0074 


Nr =N*1 

OOOOlOrO 

007^ 


1.41 TO ^ 

00001030 

0077 


1 = N ♦ 1 

0 0 001 Oho 

0074 



OooolObO 

007v» 

S 

• *i^ n r = i»iL 

onooiObO 

OOhO 


YY - f ( r ♦ 1 *U) “Y * 1 .N) 

M'lnoio/u 

OOM 


//=/ (1 ♦ i *N) - 7 ( l • ^) 

00001060 


ll 

M ( l*M =s jw| 

nncMM osu 

00^4 ^ 


I L c = n * / 

onou 1 1 00 

O0H4. 


IL J-IL* < 

oooou lu 
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00«*% 


IL4»IL*4 

OOOOlUO 

008^ 


Y (lL?*Nl«Vl?fNI 

00001130 

00H7 


VnL3%M*Y 

00001140 

00H8 


Yat4*fO*Y|'**M) 

oooomo 

00H9 



000Ullt>0 

00<>0 



OOOOllfO 

OQQl 



ooocnitiu 

00^? 


l-L'tlUP*N>=FUl«Nl 

00001 190 

0093 • 


FL<1L?*MsFLI?«N) 

00001?00 

0094 


FL<lt3*NJsFLt3*Nl 

OOOOUIO 

009S 


Itilsl 

00001??0 

009#h 


1H?«? 

0000l?30 

0097 


l«3^3 

C0001?40 

Q09P 


IFHlLL^eO.lU 00 TO 480 

000012S0 

0099 



00001?b0 

0100 


IL/=It 

00001?70 

0101 


Ir^niLC-? 

00001280 

010? 


If^?=lLt-l 

00001290 

0103 


I^lalLL 

00001300 

0104 

480 

on 8 l3m*ll.3 

OiOUlllO 

OlOK 


i>9zn) V(^n^i,N)-?ti*Nn/FC(T*N) 

00001320 

OlOr. 

H 

ij8v ( n ( y a ♦ 1 •^4) «y ( I •n)) /f l u «m) 

0000X330 

0107 


00 9 I = 

OOOUl 340 

010^ 


o/ti)5<a^iu-n*FLn*N)^'%^/f i)oFi (1-^ *Nn/(Ftn*r«>«Fcn-i«Mi 

00001380 

0109 

9 

OYU) = (nr^Y t I-inFLC l«N) ♦OMY n )*FL ( l-l •H))/ CFL U*M) ♦FL(I-l*M) ) 

OOOOliOO 

Olio 


1)0 1 4 

00001370 

0111 


04? ( n r (o? n ♦ n -o/ n n /fi. < i ^ni 

OOOOllnO 

Oil? 

13 

ti ) r CD? CP i )-ov cm /kl ? ! 

00001 ?40 

0113 


r»u l«lr3*lL? 

00001 400 

0114 


= (I-l )<>FlU»^4) ♦»>r»/c n^FUC I-l f'^in/CFLUtN) ♦FLCnWN) ) 

OOCOl4lO 

Oil*- 

14 

Y*-f n = C4iUYU•l)<>^LC^•^M)*n•^VCl)4FLCI-l•N) )/CFLU*N) ♦FLCI-1»N) ) 

00001420 

011^ 


?4cn=/(PN)*')?cn^,s»PLcpN) c*mc.?>«FLUfN))*<>? 

00001-.30 

0117 

18 

Y>*( n =Y ( nn) c 1 ) c 1 *M) c n » ( n*Ni n*? 

O000l4*»0 

Ollf- 


IFULL.K9.IU) 80 10 

C00014SO 

Oll^ 


n.L^-lLL-l 

00001 4h0 

Ol?0 


p C I4 vv|^P^K0, 1 ) 60 TO 4H6 

oooun/0 

Ol?l 1 


/«>(IlL)=0#0 

D000l4>)0 

01?? ] 


/**tlL) =0.0 

00001490 

01?3 


00 *.>»4 I = PILLM 

ooool^ou 

Ol?4 


/*-c n =-i?^>nL-i) 

oouomo 

01?^^ 


YPcn^^^uL-n 

oooon^o 

01?^ 


(.O TO 

OOOOl8'iO 

Ol?7 

4HK 

VkUlL) =0.0 

00001 8«.0 

(n^f^ 


lL)=0*O 

OOOOlbSO 

01?9 


irD I = PIILM 

000U18O0 

Olio 


YMCD.-r-^cn -U 

CM»00l >70 

om 

4f<r 

?4(n=7^UL-i) 

oonon^o 

on? 


♦40 10 •.9'f 

000U18V0 

oni 

4‘V 

COMlNtU^ 

0 0(iOl<MiO 

0134 


V.^ ( 1 ) UL4) 

OOOOptlU 

013*^ 



0 000 pi?0 

on^ 


/4(j)-/>^(iL4) 

nnno 1 f>:?u 

0) 37 


/Xi') r/*>(lL'') 

OOOO lt>HO 

01 

<»av» 

ror-f 

OOoO 1 n->U 

0134 


1)0 ii> l = l*lL 

0 pool 4^0 

Ol4n 


f y = T CPI .N) -Y c P M) 

00001^70 


44 



V .'vL FACiH IS 

ub' QUAUTY 


out 


7/«/(U1*NU/(1#N1 

OOOOUmO 

ou? 


STui«7Z/fi.n*Nl 

OOOOlhVO 

0143 


CTtl)stT/>^LCl*Nl 

uOOOl 700 

0144 


UVV#s(5) 

00001 710 

OUb 


♦*<? 

OOOUl fcO 

QU^ 


SVMM sSYf>t*^Ut 

000017IU 

0147 


S/rM*^-) aSZblM) ♦V/%»«YV 

00001740 

0|4^ 


^^•Aysi^r^rdM 

OOOOlTnO 

014V 


U(crU)*Gt.u.O) GO ro r 

OOOOl 7hO 

OISO 


t i \ *S> «Pl -4^ 

(10001770 

oisi 


0-0 ro p) 

OOOOlYtlO 

OlH? 

7 

rc r »'^UA\ 

00001 7V0 

0153 


!► ('^n n •Lf.O.O) 

oooauuo 

01S4 

10 

^o^^NUf 

OOOOlfllO 

0155 


r^ni^'UY^'cii 

OOOOU^u 

015#% 



oooouio 

0157 


TilL'"*N)s?ll#ig| 

0O00U40 

015>* 


»>0 #> I *l • ll. 

onoouso 

Olsa 


u-ui 

0000 Ut>U 

OUfi 


T^T«rn^.NUTa*Ni 

O0OOU7O 

OUl 


IF n» T. 01. ►>!►') T(rPU)sTlI^>.N)-PU 

OOOQUciO 

ou? 


ir<T'-*T,LT.U*/U) 1 Tt U*AJ)sTn^»%MUPl? 

oOOOU'^0 

0U3 

0 


OilOOUoO 

0U4 


fit ( 1 *5 > ->-t ( IL#* • 5 r *p 

OOOOUlu 

015W 


0»> ^ UUlL 

ooooivo 

Olfr^ 


n^n>a<YU*\i)-Yn*N))«$iN(^^ n*N>)*t/fi«Ni^*2ii«N))«cosrH^ n.M) 

OOOOUJO 

OUT 


S (N) 

OOOOU40 

01f\#* 



onoo 1 9 ^o 

DU^ 



uouOUau 

0170 


P‘<( n.->) 

OOOOU70 

0171 


XJI^I ,/ f < (M)-X flM^l 1 

OOOOUnO 

(i\7? 


IF(N.'^‘ . n 

OOOOU'40 

0173 


U ( XXMSl./(A(\4^l UX(N)> 

OOOO/’UUO 

0174* 


00 <f <> 1 = 1 •U. 

oooo?aio 

017^ 


(r,^^ < I ) .jv' ( !♦ P 1 

OOOOi-ar'O 

01 7*- 


l^(>^.eo•^*u lio ro u 

• MlOOf^U JO 

0177 


( I 1 

OCOO^-OhU 

\>\ 7 u 1 


ON^ 1 ! ) =0^4>AX'^ 

0000? Jt>0 

017^5 \ 


P (vi.ro. 1 ) (W* TO U 

0000?0f>0 

OUO \ 


(»%( r .M (f)N< ( I 1 • > .AM <»XX«XXP 

OO0O?O /() 

OUl ^ 

1 

t.o ro u 

u 000 ?OhU 

ou? 

\? 

PN 1 1 ..•♦) 2>>\> t \ ) 

(tnno?0'»u 

OUl 

u 

•)\ ( r *r » n ( 1 ) -Mt t<»st ( n 

OnOrtr-piO 

0 1 M 4 


HA-:. ( r ) r.'i'i ( J .M '».S 

f 1 0 0 0 ^ 11 0 

nuH 


po ?<* I=l*ll, 

OOOU? Uu 

0 \ 


Y t r o> «Y n .N) 

(lOftOr 1 

0 U 7 


/ / s / ►- ( j r - / ( r • A ) 

0)0 0?1hO 

0 Uh 


hS(j. n = 'i fvy«vf ♦//«//) 

FOOO?UO 

0 1 


As^ / r if/ //• s (.> . n ) 

oono? loO 

01‘»c 


U ( T Y , f» ,•>) 00 TO u 

o0O0?l 70 

OUl 


71 

If 

1 

p 

./» 

0000? UO 

nu<^ 


iM> r«) ir 

0 0»U)>* Uu 

OUl 

IS 

l.'A s 

0 OoOr'^uU 

0 u*« 


!► (^/.L T.u.^o 1 r’ ♦ 4 ^. 

(lOnOr ^ ! 0 

OUK 

17 

-'s ( j. r mg 


0 1 


IF . T ( 1 ..\) )M> ( J, 1 ) 

00 0 


45 



0147 

20 

OSCjtliaG 

00002240 

014A 


K<5-1L 

000022!>0 

0144 


00 JO jal.lL 

000022b0 

0200 



00002270 

020] 


wSf J*fLP)SKi><Jtll 

000U22n0 

0202 


*<S=Hb^ 1 

00002^0 

0203 


4=KN 

OOO02J00 

0204 


Ufi JO i=uic 

00002J10 

020S 


sTT5?;T(ji4CT<n-cr(j)4STin 

00002120 

0206 


KtfK^lL 

00002330 

0207 


o=i)s(j*i4n 

00002340 

020M 


If «».GT*Ta*Mn *^H536 

000023!>0 

020V 


ir •f- ) 

00002360 

0210 


<.>*L T.T (!••'*) ) f>H>sO^Pl2 

00002370 

0211 


CTT=CT<j)4cr (T) ♦sr (ji^snn 

000023HO 

021? 


os(j.n n 

0O002JV0 

0213 


ifij.j^o.n L»srj*i)=-Pi 

00002400 

0214 


i4HL-AL06(*-S(.J*l^ II/45CJ*!) 1 

00002410 

021*^ 


ssii • « j# 1 1 =c r r«M^t,-os ( j* n *STT 

00002420 

0216 

JO 

A,|I (H 1 = w»?l«STT*OS( t)'»CTr 

OOO024 jO 

0217 


C4LU (aj: • It •oo*L*"^) 

000024.«0 

02l» 


call 6**Hr-.J(K il tON'ttSiviLt lL«n 

00002430 

021'* 



U000<4bU 

0220 


<*.*-1 Tf (k; ^00) 

00002470 

0221 


(-.•S03) (Sini «Isl*ILI 

000024^10 

022? 

SO 

Ou K| jrl*lL 

OOOU24VO 

022 ^ 


y^e < J*N) =0,u 

00002600 

0224 


tj) 

00002* >10 

0?2*- 


. ,)(, I ♦ iL 

'»a0026^0 

022k 

*?v 

= 1/* < j.'O ♦St ( 1)4*. Jf n 

000026JO 

0227 

b7 

J*.a) =0\ ( J. ■) (2,442( J*l<) 1*42 . • 

OO002S4(i 

02?^ 


!► ( Ir^'-.hu. 1 ) «'0 ro ^0 

00002bsll 

022'- 


^U) s(‘ I - 1 • ! L 

00002660 

0210 


.YP^>^>=YP(^ 

Oi)002S7fi 

0211 


/* f ( I ) -r/P(D 

000o2'5rt0 

0212 


iK» J-I-IL 

0000<^^*-0 

0231 


r f j. ? 1 =!)'»( j* f) 

OOOUf'KOO 

02 U 

TO 

YPi^f ll 1 ) 

0n002e>lit 

0 236 


rtp) =^2ni) 

00002i>^0 

0236 


IF ( I >^PP.r:0#o) TO lOO 

U00U2hJ0 

02J7 


if, MO t=lflL 

00002^40 

021H 


t o *^o J*=l'. iL 

0 OOO/k^^O 

02 S'* 


// = /-«- Cv/) -/ { T .fi) 

Ot)0 02K»iU 

0240 
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PARTm 


OTRIGI'aA^ 

OF POOR OOAU 


MODIFICATIONS FOR CROSSECTIONS WITH CORNERS 
A. Dlscnssion 

Parts I and II describe a program to compute force coef^cients and 
pressure distributions over arbitrarily but smootbly shaped crossections in 
Sie absence of corners. Although solutions based on slender body theory 
are invalid over regions of high surface curvature they are still capable of 
yielding good results away from such irregularities provided additonal care 
is exercised in the computation of geometric surface properties as a corner 
is approached. Analytically, a corner represents an arbitrary break in the 
structure of local surface properties. Any scheme of specifyii^ corner 
pri^rties by a finite number of discrete parameters must involve implicit 
assumptions regarding the behavior of such corners between points at which 
data is given. For tiiis reason it is desireable to have a procedure which 
allows file user some discretion regarding fiiese assumptions without re> 
quiring an excessive amount of data to define surfaces. In the following 
procedure this discretion is excersised in the choice of the distribution of 
ortiiogonal lines S| introduced in Fig. 2. 

In a finite computational scheme the difficulties inherent at a cormr 

first become manifest when the local radios of curvature on C becomes 

n 

small compared to file distance in fiie y, z plane to the neighboring cros- 
section Such points are illustrated in Fig. 16 at (i,n) = (15,5), (15,6), 

(15,7). For practical computations such points are equivalent to the sharp 
corners of (4, 2), (4, 3) and must be treated in the same way. In contrast 
to the procedure of Part n which "rounds off" regions of higher curvature 
it is more appropriate now to adapt the opposite procedure namely: a region 
of finite but large curvature is to be rejdaced by a sharp corner. If this 
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approacimation should prov« too erode it would dfesn be necessary to include 
an additional contour between C and C . , so that tfie distance between 
contours is less fliaa tbs local radius of curvature, a procedure uddeh is 
equivalent to supplying more detailed date to fill in tee objectionable gaps. 

B« Definitions 

1. Stringers S| 

Tlie lines orteogonal to the croesecfional contours C(n) and for wfaicb 
i s constent shall be called stringers. These are tee family of lines first 
illustrated in Fig. 2. 

2. Corner lines i s IC(K) 

These are lines passing t h r u corner paints. They are to be considered 
as part of tee family of strii^ers S.. As such they are continued over tee 
entire lengte of tee body even though previous or .subsequent crossections do 
not have corners. An example of one such line is shown for i s 4 in Fig. 16. 

Corner lines are distinguished by the index IC(K) s i signifying teat the index 

fh 

of tee K corner, counting counter clockwise, is i. Thus in Fig. 16 
IC(2) = 4. (For programming convenience it is expedient to designate the 
first stringer is 1 as a corner line ie (C(l) = 1 even though teere may be no 
corners along this line. ) 

3. Submerged lines IBP(K, n). IBM(K, n) 

A stringer S| from tee contour C^ may intersect a corner line before 
it intersects tee next contour as illustrated in Fig. 16 at (10,6), (16,4), 

(14,5) and (17,6). Subsequently such stringers are considered to follow the 
corner line and are regarded as submerged. At the corner on C^ the 
highest submerged line index is denoted by IBP(K, n) and tee lowest by 
IBM(K,n). Thus from Fig. 16 we find IB1(5,7) s 17, IBM(5,7) = 14. A 
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coriwr lias may also ba coaoted as a submerged line ie: 1BM(15, 5) s IS, 
1BM(15,4) s IS. We note ttnn, tiiat every intersection of a corner line 1C(K) 
wlUi a contonr bas associated widi it tiie indices 1BM(K, n), IBP(K, n). 

For porposes of illnstratlon a cc...plete table of 1BP(K, n) is provided in 
Fig. 16. Finally, we note tiiat in tibe absence of any corners along i s 1 we 
net lBld(l,n) » XL -f 1. In Fig. 16 this means Biat IBM(l,n) s 19 . 

C. Modificationb to flie compotational procednre 

1. Collocation Points 

Points P'(i,n) at which dv/dx, a etc. are to be evaluated were pre- 
viot«Hly found by smoodi interpolation between P(i-2,n) and P(i-l-2,n). To 
avoid dw requiring of an excessive number of data locations P(i, n) between 
corners this has been modified so that P '(i, n) is read directly from supplied 
data or by simple interpolation between Mighboring locations P(i,n), P(it-l,n). 
In many practical applications the contour curvature between two corners is 
small and die later procedure should be adequate. 

2. Computation of 5y/^ 

Values of 6 v/ 6 x are to be found at P(i, n), P(i-l-l,n) and interpolated to 
obtain a value of P'(i, n) between i and i -f 1. (This represents a minor but 
necessary change from the procedure of Part II which determines 6 v at 
P'(i, n-1) 8 t P'(i,n) and interpolales die associated derivatives along the x 
direcdon). The increments 6 v are taken along the stringers and as long as 
these do not intersect the corner lines the determination of 6 v/ 6 x at the data 
points P(i, n) is carried out as thc^gh no corners were present. 

When a stringer intersects a corner line the local corner geometry 
is assumed as shown in Fig. 17 which represents an enlargement of the local 
configuration as it appears in Fig. 16 at P(4, 2) and P(4, 3). While &v as 
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indicatad in Fig. 17 may b« calctdated directly from tiie data presented in 
dm plots of and S|, flie value of be most be inferred from die assump- 
tion diat die corner line shown in Fig, 17 is closely approsimated by a 
straight line, Tfans witii 6 v 2 » as indicated in Fig. 17: 

This is to be compared widi the calculadon away from a corner wl»re 6 (x) 
is simply x(nl-l) - x(n). 

To devise a program which is applicable to all possible instances of 
corner geometry it is necessary to have tests which indicate when a stringer 
emerges from corner as between P(4, 2) and P(4, 3) in Fig. 16, and when it 
converges toward a corner to become snbsequendy submerges as is the case 
between P(ll, 6 ) and P(ll,7). Such a test is readily constructed with the aid 
of the indices IBP and IBM. Thus for example: 

at least one stringer 

IBP(K,n-fl) < IBP(K.n) has emerged between 

P(IC(K),n) & P(IC(K), n+1) 

and 

1BP(K, n) - IBP(K,n-fl) = no. of emerged stringers. 

In this manner IBP and IBM provide complete information regarding the 
emergence or convergence of stringers on either side of a corner line. 

This information together with implied geometry of Fig. 17 enables the 
computation of bv/tx at the center of contour segments which are adjacv.nt 
to cori^r lines. 

3, Curvature 

The fact that curvature is divergent near corner -like points leads to 
errors in the computation of dep/dx when using the program of Part II. This 
program in effect rounds off corners whereas as pointed out in the discussion 
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abovtt a mora appropriate procedure is to treat regions of high curvature as 

sharp corners ie as though high curvature regions were concentrated at a 

corner point. Witih this procedure the curvature of segments adjacent to a 

corner is small and may be obtained by extrapblatton from a neighboring 

segment. Thus, referring to Fig. 16 we would have: 

h(3,3)sh(2,3) 
h(4.3) = h(5.3) 


4. Computation of 6q/oat 

In Part U 6o/6x was approximated by central divided differences 
involving o(i,n-l), o(i,n), o(i, n-fl). This procedure breaks down at a cor- 
ner. The rules to be followed near corners will now be that (6(t/6x). will 

1 , n 

be computed by: 

Forward divided difference when S. and/or emerge from 

a corner. 

Backward divided differences when and/or converge 

to a corner. 

Central divided differences away from cc *ner. 

As an illustration corresponding to Fig. 16 

In the evrnt o£ a stringer emerging just behind a segment and again 

converging just ahead of it ^G/dx shall be assxuned to be zero. 


5. dG/ds 

To compute dG/ds we just find G at all the data points P(i,n) (except 
at a corner pi. ) by interpolation between neighboring collocation points 
P*(i-l,n), P^(i, n). At corner points dG/ds is then found by forward 


57 



differences leaving a corner along and by backward differences when 
apinroaching a corner. 

6. Matriac Inversion and Summation 

Tor the matrix inversion process encountered in the evaluation of 
a(i,n) it ^s convenient to reorder tiie indices so that the actual finite segments 
of a contour are indexed consecutively. This involves shipping over sub« 
merged segments in the counting process. Such a reordering may be accom- 
plished fiirough file introduction of a new index IR(m) for which the m are 
consecutive indices and: 

lR(m) s i 

for values of i correspondii^ to unsubmerged segments. Thus we would have, 
for example 

mL 

^ c(lR(m), n) a(lR(m), j n) s ^ff(i* n) a(i, j, n) 
m=l 

where the latter summation is taken only over fiiose values of i corresponding 
to segments which are not submerged. 

The remaining computational procedures from Part II are not affected 
by the presence of corners and do not require modification. 
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FIG. 2 CROSSECTION BOUNDARY SEGMENTING 

SCHEME IN BODY AXES COORDINATES 
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FIG. 3 DETAILS OF VARIABLES PERTAINING TO 
SEGMENT i,i+l OF BOUNDARY Cn 
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FIG. 4 RELATIVE POSITIONS OF Cn AND Cn^l 
IN BODY AXIS AND WIND AXIS 
REFERENCE FRAMES 




FIG. 6 «^/4x AT X -I ON DRCULAR CONE AT 

ANGLE OF ATTACK.a>«.l 











FIG. 8 LIFT COEFFICIENT FOR QRCULAR 
CONE AND OGIVE AT ANGLE OF 
ATTACK, a -0.1 
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F*'-. 9 MOMENT COEFFICIENT FOR CIRCULAR 
CONE AND OGIVE AT ANGLE OF 
ATTACK, a -0.1 
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FI6.K) SAMPLE FUSELAGE 


68 








FIG. 12 Cp ALONG UPPER AND LOWER 

MERIOAN OF SAMPLE FUSELAGE 
AT ANGLE OF ATTACK, o»IO» 
AND MACH NO. >0 
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FIG. 16 ILLUSTRATION OF SEGMENTING SCHEME 
FOR CONTOURS WITH CORNERS 
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